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Abstract 

A system of two coupled cavities is studied in the context of bipar- 
tite, continuous variable entanglement. One of the cavities is pumped 
by an external classical source that is coupled quadratically, to the 
cavity field. Dynamics of entanglement, quantified by covariance mea- 
sure [Dodonov et al, Phys. Lett A 296, (2002) 73], in the presence 
of cavity-cavity coupling and external pumping is investigated. The 
importance of tailoring the coupling between the cavities is brought 
out by studying the effects of pump fluctuations on the entanglement. 
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1 Introduction 

Quantum information processing requires entangled states [TJ 12] ■ Many 
physical systems such as Josephson devices, trapped ions, NMR, have been 
considered for quantum information processing. All these system are en- 
dowed with finite number of states for manipulation. Of late, systems which 
are equipped with infinite dimensional Hilbert spaces have been analyzed 
in the context of teleportation, computation, etc [5J HI [SJ El E] • Nonlinear 
optical processes, cavities in photonic bandgap materials are capable of pro- 
ducing states of light that are useful in quantum information processing. 
A notable feature of bosonic modes is their robustness in retaining their 
entanglement in spite of thermal noise [8]. Also, experiments with entan- 
gled photon states are of significance in testing the foundational aspects of 
quantum theory [9j. 

Microcavities have been extensively employed in generating atomic qubits, 
engineering of quantum states, creating polarization entangled photonic 
states to do proof-of-principle tests such as verification of Bell-type in- 
equalities, entangling atoms and/or photons [13 [HI H] , studying the phe- 
nomenon of decoherence|13j . etc. More interestingly, coupled cavities have 
been thought of as a conduit to transfer and engineer entanglement between 
distant atomic qubits [H], which would be required for distributed quan- 
tum computing. Cavities containing atoms (either positioned or on-flight) 
have been used in the previously mentioned contexts. It is possible to use 
suitably prepared atoms to modify the cavity field so as to make the later 
useful in applications such as photonic qubits [TU [16], optical Josephson 
interferometers [T7] . etc. 

In this work, a system of two coupled single-mode cavities is considered. 
The cavities are coupled by evanescent fields. Hence, the coupling strength 
can be tailored by the proximity of the cavities and the permittivity of the 
cavity walls. The field modes of the two cavities form a bipartite, bosonic, 
continuous variable (CV) system. Because of the coupling, the modes of the 
two cavities can get entangled. In addition to the inter-cavity coupling, one 
of the cavities is externally pumped. The results presented here pertain to 
the case of treating the pumping field classically. 

The Hamiltonian of the system under consideration is 



The operators a and a* are respectively the annihilation- and creation- op- 
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erators of the mode (a-mode, hereafter) in the cavity which is externally 
pumped. Similarly, b and b' are the relevant operators for the mode (6- 
mode, hereafter) in the other cavity. The strength of the coupling between 
the cavities is characterized by the parameter A. The magnitude of this 
parameter is related to the hopping strength of photons to hop between 
the cavities. The term ea) 2 + e*a? corresponds to the contribution to the 
Hamiltonian from the interaction between the external pumping and the a- 
mode. This interaction, for instance, could be achieved by an active medium 
within the cavity. The magnitude of this interaction strength between the 
a-mode and the classical driving field is given by the parameter e which 
incorporates the effects of the active medium and the external pumping. 
A similar approach is adopted to explain the maser action which includes 
many external modes [H]. The external classical field affects entanglement 
as it delivers energy into the pumped cavity which supports the a-mode. 
However, whether the entanglement is enhanced or decreased depends on 
the coupling (A) between the cavities. Firstly, the dynamics in the absence 
of the external pumping is studied. These results are useful to understand 
the effects of pumping. 



2 Dynamics of fields: without pumping 

The Hamiltonian to study the dynamics in the absence of external driving 
corresponds to setting e = in H. In that case, the Hamiltonian evo- 
lution is easily obtained since the commutator [a^a + b'b, a)b + aW] van- 
ishes. This implies that the evolution operator exp(-itH) can be factored 
as exp(— it(a)b + ab')) exp(— it{a)a + »'&))■ This type of Hamiltonian is 
known in other contexts such as the symmetric beam splitter [19], the para- 
metric conversion in a medium with oscillatory dielectric constant [20J and 
has recently been used to examine intrinsic phase coherence in a laser [21]. 
The bilinear coupling \(a^b+ab'), between the modes cannot entangle them 
during evolution if both the modes are in classical states initially [22| 123]; for 
instance, the a-mode in a coherent state and the 6-mode in vacuum state. 
To generate an entangled state from a product state, during evolution in 
the presence of bilinear coupling, at least one of the modes has to be in a 
nonclassical state. So, it is meaningful to assume the initial state of the 
cavities to be a number state. Since all number states, except the vacuum, 
are nonclassical, the bilinear coupling may entangle the a- and b- fields. It 
is experimentally possible to prepare the cavities such that a-field is a Fock 
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state, say \N), and the 6-field to be vacuum state. 

Another consequence of the vanishing commutator [a'a+b'b, aJb+alr] = 
0, is the existence of nontrivial, irreducible, invariant subspaces. Each in- 
variant subspace is the span of the states of the form {\m,N — m)}, where 
m runs from zero to N. The state \m, N — m) of the modes in the cavities 
represents the situation in which the a-mode is in the number state \m) and 
the 6-mode is in \N — m). There is an invariant subspace associated with 
each value of N, an eigenvalue of the total number operator a^a + b'b. All 
the states in an invariant subspace are the eigenstates of the total num- 
ber operator. Further, the invariant subspaces corresponding to different 
total quantum numbers are disjoint. Taking the initial state as |iV, 0), the 
time-evolved state \ip(t)} = ex.p(—iHt)\N,0) is 

N 

\ip(t)) = exp(-iNujt) cos N (Xt) 

n=Q 

This state belongs to the invariant subspace corresponding to a total quan- 
tum number N which contains the initial state | A, 0). An interesting feature 
is that the states of the form \m, N — m) are the only product states in the 
invariant subspace. All the other states in a given invariant subspace are 
entangled [24J. The coefficients in the Fock state expansion are the bino- 
mial coefficients and the states defined in Eq. [2] are referred as two-mode 
binomial states [25]. 

In this work, the covariance criterion is adopted as the measure of 
entanglement [26] . It is expressed as 

|aftt| 2 + \hb\ 2 
[2(^+1) (Sts + |) 

The bar is used to indicate covariance, for example, ab> = (ab') — (a)(b'), 
where (• • •) stands for quantum mechanical expectation value of the relevant 
operator. The quantity Y is non-negative and less than unity. For product 
states, the value of Y is zero. Nonzero values of Y implies the state is 
entangled. But there exist entangled states for which Y is zero and hence 
the criterion is not universal. Nevertheless, it is easy to compute and useful 
in identifying and quantifying entanglement when it assumes nonzero values. 
If the initial state is of the form | N, 0) , the covariance parameter of the 
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evolved state given in Eq. [2] is 

Y = N\sin(2Xt)\ ^ 

2 [2(iVcos 2 (At) + .5) (TV sin 2 (Ai) + .5)] 1/2 ' 

With nonzero coupling, the maximum value that Y attains during evolution 
is N/\^2(N + 1) and it is attained whenever At is an odd integral multiple of 
it/ 4. It is clear that the peak value of Y attained during evolution increases 
with N and approaches l/v2 asymptotically as N becomes large. Further, 
the peak value depends on the total quantum number N and not on the field- 
field coupling constant A. In Fig. Q] the variation of Y with time is given 
for different values of the total quantum number, which is the eigenvalue of 
the operator a) a + b*b. The time of evolution is measured in the units of 
7r/A, which is the temporal periodicity of pumping- free evolution. In all the 
subsequent discussions, whether the evolution is periodic or not, the time is 
always expressed in the units of tt/X. 

The universal measure of entanglement for bipartite pure states [27] is 
the von Neumann entropy S defined as 

S = -TV [p a log 2 p a ] , (5) 

where Tr stands for trace and p a is the reduced density matrix for the a- 
mode field. It is instructive to study the behaviour of von Neumann entropy 
(S) vis-a-vis Y in the present case. For the states defined in Eq. [21 the 
reduced density matrix p a for the a-field is 

N W 

p a = cos 2N (Xt) £ ' tan 2 «(At)|iV -n)(N-n\. (6) 

^— 1 nl (JS — n)\ 

Denoting the coefficient of \N — n)(N — n\ in p a by c n , the expression for the 
entropy is S = —2J2n=o l c n| 2 log2 \ c n \- In Fig. [2j the variation of S and Y 
are shown as functions of time taking the initial state to be |5,0). The two 
measures exhibit very similar features. When S is nonzero, the covariance 
measure Y is nonzero too, indicating that it is a good criterion for the kind 
of entangled states given in Eq. [2j The peak in the entropy occurs when 
the two modes have nearly equal number of photons, that is, (a^a) ~ (b'b). 
Hence, entropy peaks correspond to minima in the difference of the photon 
numbers of the two modes. The photon number in the a-mode is N cos 2 (Xt) 
and in the 6-mode is ./V sin 2 (Ai). Therefore, (a^a) equals (b'b) whenever At is 
an odd multiple of 7r/4. The reason for the increase in the entanglement as 
the total number of photons increases is readily inferred from the expression 
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for entropy S given after EqJBJ The entropy S is maximum, subject to the 
constraint that Trp a = 1, when all \c n \ are equal. There are N + 1 terms 
in the expression for S and the maximum S is log 2 (A r + 1) which increases 
with N. 



3 Effect of external pumping 

The evanescent coupling allows photons to hop from one cavity to the other. 
This allows a redistribution of energy between the cavities. Since the o-mode 
is subjected to external pumping, the energy of the pump field is fed to the 
a-mode which, in turn, is used to energize the 6-field. The pump field is 
treated classically and hence the driving term in the Hamiltonian is a func- 
tion of the operators of the a-mode. For weak couplings, it is reasonable 
to assume the function to be linear in a) and a. Classically, this amounts 
to an interaction proportional to the product of the field strength of the a- 
mode and the amplitude of the driving field. If the coupling is stronger, 
the function must involve higher order terms involving the creation and the 
annihilation operators. In what follows, the effects of quadratic coupling on 
the entanglement between the cavities are discussed. As a remark, it is noted 
that if the pumping is linear, that is, the driving term is ea) + e*a, there 
is no effect on the value of the covariance parameter Y. There is indeed a 
change in the dynamics and, consequently the entanglement is affected. But 
the dynamics of the covariance criterion is insensitive to such modifications, 
a reminiscence of the fact that the covariance criterion is not universal. The 
quadratic coupling, however, leads to significant changes in the dynamics of 
Y. 



The Heisenberg equations of motion for the operators of the two mode 
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The equations imply that the evolution of the operators of the two modes 
are coupled. The coefficient matrix has no explicit time-dependence and 
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hence the equations are easily solved to obtain 
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where M represents the coefficient matrix in Eq. [7J 

The operator exp [—itM] can be expressed in terms of lower powers of 
M using Cayley-Hamilton theorem[28|. Let the four eigenvalues of M be de- 
noted by a, f3, 7 and 5. Defining A = cu 2 +X 2 -2e 2 and B = \/lu 2 \ 2 - X 2 e 2 + e 4 
the eigenvalues are 

-P = \J A - 2B, (9) 
-S = VA + 2B. (10) 
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These are the eigenvalues if e is real and they are distinct if the determinant 
of M is nonzero, which is always the case for the typical values of the 
parameters assumed. Now, applying Cayley-Hamilton theorem, 

'sin at sin7t s 



4Bexp 



-itM 



7 



M A 



+ (cos "ft — cos at) M 2 



a sin 7t 7 sin at 



7 



Q: 



— ^7 cos at — a cos 7* ) I4, 
where I4 is the 4x4 identity matrix. 



M 



On using the matrix exponential, refer Eq. [TT| in the solution given 
in Eq. [8j the expectation values of time-evolved operators are expressed in 
terms of their initial expectation values. Let the time- dependent coefficients 
of M 3 ,M 2 ,M and I4 in the expression for exp[— itM] given in Eq. [TTj be 
C 3 ,C 2 ,C\ and Co respectively. Further, define 

U U) x = C + wCi + (A - 2e 2 )C 2 + w(2A 2 - 2e 2 + A)C 3 , 

K,,A = XC 1 + 2XuC 2 + A(2w 2 -2e 2 + A)C 3 , 

W^x = 2eCi + 2e(A 2 + A)C 3 , 

X^x = Co + cod + (A 2 +lo 2 )C 2 + cj(3X 2 + uo 2 )C 3 , 

Y^ x = 2e(XC 2 + coC 3 ), 

Zuj,x = -2X 2 eC 3 . 
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If the initial state of the two cavities is |iV, 0), the relevant quantities to 
compute Y are 



aSt = = {l + N)U^ x V^ x + NW^ x Y_^ x + V^ x X_^ x , (13) 

ata = =(i + n)u^ x u^x-nw^ x + v^ x v^ x -i, (u) 

Wb = =(l + N)V U] , x V^ x + NY UJiX Y^ x + X u}iX X^ x -l. (15) 



Using the expressions in Eqs. [T2lll5l for the various bilinear combinations 
of the creation- and annihilation- operators of the two cavity modes, the co- 
variance criterion of entanglement is computed. In Fig. [3] the evolution of 

Y with time is shown for various combinations of A and e. When both A 
and e are chosen to be 0.1, the entanglement measure Y builds up to 0.6. 
The occurrence of peak values in Y is corroborated with the occurrence of 
minimum in the difference between the photon numbers of the two modes. 
The variation of the ratio of the mean photon number difference \{a^a — b'b}\ 
to the total number of photons \(a*a + Wb) | is shown in Fig. [H If the ratio 
is unity, the photon number of one of the modes is zero. If the ratio is zero, 
the two modes have equal number of photons. The evolution of the ratio 
is shown for three different combinations of A abd e. If the cavity-cavity 
coupling A is sufficiently strong that the energy increase due to the external 
pumping gets distributed to both the cavities, the difference in the average 
photon numbers in the two cavities can become zero, as in the case of A 
and e both being 0.1. If A is small, the energy transfer between the cavities 
is ineffective and the photon number difference is large, for instance, when 
A = 0.001 and e = 0.1. To illustrate the effect of larger A, the evolution 
of the ratio if A = 0.005 and e = 0.1 is shown in Fig. |U In this case, the 
difference becomes smaller than the e = 0.001 case; however, the coupling is 
still not strong enough for efficient energy exchange between the cavities to 
make the difference in the photon number to become zero. 

The occurrences of vanishing photon-number difference and the peak in 

Y happen at same instants, seen by comparing Figs. [3] and [H When the 
pump-field coupling is large, the o-mode gets energized at a faster rate. If 
the coupling A is such that the exchange of energy between the cavities is 
effective, the a-field and the b- field attain nearly equal energies so that the 
difference of their photon numbers becomes small at various instants. This, 
as in the case of pump-free evolution, is correlated with the occurrence of 
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(12) 



S 



entanglement peaks. Similarly, if both A and e are equal but smaller, chosen 
to be 0.001 for the purpose of discussion, the slow rate of energy build up 
in the a-mode is matched by the flow of energy between the modes. Hence, 
in this case too the energies of the two modes can be nearly equal during 
evolution. At such instants, the entanglement becomes higher. 

If (A << e), for instance, A = 0.001 and e = 0.1, the external pumping 
primarily enhances the energy of the o-field as there is no effective flow of 
energy to the 6-field because of low value of the coupling between the cavi- 
ties. This increases the energy difference between the modes. The evolution 
of Y in this case is shown in Fig. [3j Compared to the pump-free evolu- 
tion of the initial state |5,0) shown in Fig. [H the maximum entanglement 
attained is smaller as the energy difference between the modes is larger. 
The 6-mode is initially in vacuum state. Due to the weak coupling between 
the cavities, the 6-mode remains in a state in which only the vacuum state 
and a few number states are present. Therefore, the entanglement becomes 
smaller. The entanglement attains a peak whenever the energy difference 
is minimum. Thus, generation of large entanglement requires that energy 
flows into the 6-field, which, in turn, requires a strong interaction between 
the cavities. 

In the opposite limit A >> e, for instance A = 0.1 and e = 0.001, the 
total energy of the system does not increase much as pumping is weak. The 
dynamics is very similar to that in the case when there is no pumping. 
The effective coupling between the modes leads to nearly equal energies in 
the modes during evolution. Whenever the mean photon numbers of the 
modes match, the entanglement attains a peak. The evolution profile of Y 
is not distinctly seen in Fig. [3] as it merges with the profile corresponding 
to A = 0.001 and e = 0.001. 



4 Effect of pump fluctuations 

The maximum value that Y attains during evolution is dependent on the 
initial photon number N, the coupling A and the pump-field interaction 
strength e. In Fig. [5] the maximum value of Y attained during evolution 
is shown as a function of e. Results corresponding to two different values 
of A are given taking the initial state to be |5,0). It is seen that the peak 
value increases with e if A is 0.1 whereas it decreases drastically if A = 0.001. 
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If A = 0.1, the peak value of Y is close to 0.6, for small values of e. This 
is nearly the value that is attained in the absence of pumping since the 
pumping is weak. As the value of e increases beyond 0.05, the peak value 
increases. For sufficiently large e, say, « 0.4, the peak value corresponds 
to the maximum attainable in pumping-free case and initial N as large as 
50. This is understandable as large pumping leads to large energy input. 
Since A = 0.1 is sufficient enough for energy transfer between the modes, 
both the modes will have more number of Fock states in their superposition. 
Hence, a large entanglement is possible. If A = 0.001, a representative value 
for small A, the coupling is not good enough for the energy to flow from 
the a-mode. Thus, the fo-mode will have fewer number of Fock states in 
the superposition. Also, as e increases, the difference in the photon number 
of the modes will become larger. So, there is reduction in the maximum 
entanglement between the modes as compared with the values in the absence 
of pumping. The dependence of peak Y on e is shown in Fig. [5] for two other 
values of A. It is seen that as A becomes larger, the peak Y during evolution 
is nearly constant over a larger range of e. 

From Fig. [5J it is seen that if A is 0.1, the maximum attainable Y varies 
from 0.6 to 0.7 as e increases from to 0.5. On the other hand, if A = 0.001 
the maximum attainable Y changes from 0.6 to 0.03 as e varies from to 
0.5. The change in maximum Y is larger for the smaller A implying a larger 
sensitivity to changes in e. This has important consequences. In any pump- 
ing scheme there are fluctuations in the driving amplitude. The effect of 
these fluctuations on the entanglement is decided by the sensitivity of the 
system. From the discussion it is expected that the dynamics should not 
differ much if A is large whereas perceptible changes are possible if A is small. 
To account for the driving field fluctuations, the coefficient e is assumed to 
fluctuate about its mean value. The fluctuation is assumed to be Gaus- 
sian with mean equal to the value of assumed value of e and variance to be 
one-tenth of the mean. When the mean value of e is 0.001, the fluctuations 
do not affect the entanglement. But at higher mean values, fluctuations 
in e affect entanglement dynamics. Taking mean of e to be 0.3, the time- 
evolution of Y is shown in Fig. [6] for A = 0.001 and 0.1. The fluctuating 
driving term is taken to be piecewise constant function. The total evolution 
time is divided into 100 equal parts; in each part, the value of e is chosen 
at random from the Gaussian distribution defined earlier. The evolution is 
tracked for ten different random trials in each case. The evolution of Y cor- 
responding to small A, shown in (a) of Fig. [6j shows a visible spread in the 
evolution during different trials. In order to distinctly bring out the effects 
of fluctuations, the dynamics is continued for five units of scaled time. The 
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entanglement itself is very small as the field-field coupling is small. But the 
fluctuations induce large relative changes in the entanglement. Thus, strong 
pumping and weak coupling between the cavities lead to large uncertainties 
in entanglement. With larger A, refer (b) in Fig. [6l the spread is almost 
insignificant. Strong coupling between the modes reduces the effect of the 
driving-field fluctuations on the entanglement. As stressed earlier, sufficient 
coupling between the cavities allows for efficient flow of energy and the max- 
imum value of Y during evolution is also large. It is, therefore, possible to 
choose the cavity-cavity coupling so that the effects of fluctuations are min- 
imized. In other words, smaller coupling between the cavities necessitates a 
more stable pumping. 

5 Summary 

Two coupled cavities generate entangled, bipartite, non-Gaussian states. In 
the absence of external pumping, the cavities exchange energy periodically. 
The entanglement attained during evolution increases with the number of 
photons present in the cavities initially. Inclusion of external pumping in 
one of the cavities amplifies the field in that cavity, which, in turn, energizes, 
via the evanescent coupling, the mode in the other cavity. This affects the 
entanglement between the modes. The covariance criterion of entanglement 
is insensitive to the driving amplitude if the driving is linear. However, 
with quadratic pumping, the states occurring during the evolution are such 
that the covariance criterion identifies entanglement. If the couping between 
cavities is weak, increasing the driving field amplitude decreases the peak 
entanglement attainable during evolution compared to the pump-free evolu- 
tion. Therefore, it is important to tailor the pumping strength to match the 
cavity-cavity coupling. During evolution, maximum entanglement between 
the fields in the cavities occurs when the mean number of photons in the 
two cavities are equal. 

For weak coupling between the cavities, the maximum attainable entan- 
glement decreases rapidly as the pumping strength is increased. Though the 
entanglement increases for weak pumping, it begins to decrease for further 
increase of pumping strength. This is because of preferential increase in the 
energy of the pumped cavity leading to large difference in the mean number 
of photons of the two cavities. On the other hand, if the coupling between 
the cavities is stronger, the maximum entanglement attained during evolu- 
tion increases for larger range of values of the pumping strength. Also, the 
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change of entanglement with pumping strength is not as rapid as the change 
if the coupling between the cavities is weak. Consequently, the evolution of 
entanglement is less sensitive to fluctuations in the pumping if the cavities 
are coupled effectively. 
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Figure Captions 



Fig.l Variation of Y with time. The time-axis represents scaled-time defined 
as the ratio of time to tt/X. One unit of scaled-time corresponds to ir/\. The 
curves shown correspond to different total quantum numbers: N=l (dash), 
5 (dot-dash), 10 (continuous) and 50 (dot). In all cases the coupling A is 0.1. 

Fig. 2 Variation of covariance measure Y (continuous) and von Neumann 
entropy S (dash) with scaled time. The initial state is |5, 0) and A = 0.1. 

Fig. 3 Comparison of evolution of Y for weak and strong couplings when 
pumping is quadratic. Different plot types are used for different combina- 
tions of A and e: dashed - A = 0.001 and e = 0.1; continuous - 0.001,0.001, 
beaded - 0.1,0.1. The continuous curve hides the curve corresponding to 
the case: A = 0.1 and e = 0.001. The initial state is |5, 0). 

Fig. 4 Ratio of \(tfa - S+S> | to {a) a + 5*5) as a function of time. The 
curves shown correspond to A = 0.001 and e = 0.1 (dash); 0.001, 0.001 (con- 
tinuous), 0.1,0.005 (dots). 

Fig. 5 Maximum value of Y as a function of e. For weak coupling A = 0.001 
(continuous) and for strong coupling A = 0.05 (dots). Additional curves cor- 
respond to A = 0.005 (dot-dash) and 0.01 (dash). The initial state is |5,0). 

Fig. 6 Evolution of Y in time under fluctuating driving field. The mean 
driving amplitude is 0.3. Two values of field-field couplings are assumed: 
A = 0.001 [figure (a)] and A = 0.05 [figure(b)]. Different evolution profiles 
correspond to different trials. In each of the figures, results of ten trial 
evolutions are given. The initial state is |5,0) for all the profiles. 
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